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We consider the dynamical correlation functions of the quantum Nonlinear Schrodinger equation. 
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1 Introduction 


We consider exactly solvable models of statistical mechanics in one space and one time dimension. 
The Quantum Inverse Scattering Method and the Algebraic Bethe Ansatz are effective methods for 
a description of the spectrum of these models. Our aim is the evaluation of correlation functions of 
exactly solvable models. Our approach is based on the determinant representation for correlation 
functions. It consists of a few steps: first the correlation function is represented as a vacuum 
expectation-value of an operator-valued Fredholm determinant, second the Fredholm determinant 
is described by a classical completely integrable equation, third the classical completely integrable 
equation is solved by means of the Riemann-Hilbert problem, fourth, the vacuum expectation-value 
of the asymptotics for the operator-valued Fredholm determinant is calculated. This permits us to 
evaluate the long distance and large time asymptotics of the correlation function. Our method is 
described in |^. 

The quantum Nonlinear Schrodinger equation can be described in terms of the canonical Bose 
fields {x,t G R) with the standard commutation relations 

['ip{x,t),Tp^{y,t)] = 6{x - y), ['il;{x,t),ip{y,t)] = [ip\x,t),il;^y,t)] = 0. (1.1) 

The Hamiltonian of the model is 

R = y dx (j)x'tp\x)dx^ix) + c'tp\x)'ip\x)^{x)^{x) — h'ijj\x)'ijj{x)^ . (1.2) 

Here 0 < c < oo is the coupling constant and /i > 0 is the chemical potential. The Hamiltonian H 
acts in the Fock space with the vacuum vector |0). The vacuum vector |0) is characterized by the 
relation: 


'!/’(x,t)|0) = 0. (1.3) 

The vacuum vector (0| is characterized by the relation: 

(0|V’^(x,t) = 0, (0|0) = 1. (1.4) 

The spectrum of the model was first described by E. H. Lieb and W. Liniger |Q|, Q. The Lax 
representation for the corresponding classical equation of motion 

+ (1-5) 

ot dx^ 

was found by V. E. Zakharov and A. B. Shabat |^]. The Quantum Inverse Scattering Method for 
the model was formulated by L. D. Faddeev and E. K. Sklyanin |j^. 
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In this paper we shall consider the thermodynamics of the model. The partition function Z and 
free energy F are defined by 


Z = tr 



__F 

= e T. 


( 1 . 6 ) 


The free energy F can be expressed in terms of the Yang-Yang equation 

,(A) = ,2 + 1” (l + (1.7) 

T f e(^i) \ 

F = J Infl + e r j d/i. (1.8) 


The correlation function, which we shall study in this paper, is defined by 




tr T 

tr e T 


(1.9) 


In a previous paper |^] we obtained the determinant representation for this correlation function. 
It was shown in the paper |Q] that the Fredholm determinant thus obtained can be expressed in 
terms of solutions of an non-Abelian nonlinear Schrodinger equation. More precisely, the second 
logarithmic derivatives of the Fredholm determinant with respect to distance and time are densities 
of the conservation laws of this equation. However, such a description of the determinant is not 
complete, since the non-Abelian Nonlinear Schrodinger equation has an infinite set of solutions, and 
it is not clear a priori, which of them describes the correlation function. 

On the contrary, the Riemann-Hilbert problem associated with the Fredholm determinant is 
uniquely solvable. This is the main advantage of the approach based the application of the Riemann- 
Hilbert problem. 

The plan of this paper is the following. In section 2 we review the determinant representation 
and recall definitions and notations used in In section 3 we formulate the Riemann-Hilbert 

problem. We prove the equivalence of the integral equations considered in and the Riemann- 
Hilbert problem. The Lax representation of the non-Abelian Nonlinear Schrodinger equation is 
considered in section 4. Section 5 is devoted to the modified formulation of the Riemann-Hilbert 
problem, which is especially useful for asymptotic analysis. 


2 Determinant representation for the correlation function 

In this section we summarize the formulations and results obtained in the previous papers |Q, |Q] 
for the reader’s convenience. Our starting point is the determinant representation for the dynamical 
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correlation function of the local fields obtained in 

,-iht det^I + y^ 




bi 2 {u, v)dudv\0). 


( 2 . 1 ) 


27r det(/-^i^T)- -- 

In order to describe explicitly the r.h.s. of (2.1) we present here the system of definitions and 


( 2 . 2 ) 


notations, used in 

2.1 Vectors and operators of the Hilbert space 

Consider ket-vector \E^{X)) 

V E?(\W) 

belonging to a Hilbert space Ti.. Here E^{X\u) are two-variable functions. It is convenient to present 
the space Ti as the tensor product of two spaces H = H ^ H. Each of the discrete components 
E^{X\u), being a function of the variable A, belongs to the space TC. In turn, the space H consists 
of two-component functions of the variable u. 

In order to define scalar products in the space fi we introduce the bra-vector {E^{X)\ 

{E^{X)\ = [Ef{X\u),E^{X\u)), (2.3) 

which also can be treated as vector of TC. Then the scalar product in Ti. is given by standard formula 

CXD 

{E^{X)\E^{^i))= J {Ef{X\u)E^{fi\u)+E^{X\u)E^{^i\u)) dn. (2.4) 

— OO 

Below we shall consider two types of operators acting in Ti. The first type of operators of which 
/ H in (|2.I|) is an example, acts in the space Ti only. By definition 


/ OO \ 

Ef(A|u)+ [ H(A,/i)Ef(/i|u)d/i 


[i + v)o\E^{^i)) = 


E^{X\u)+ I V{X,fi)E^{fi\u)dfi 


Rt 


V 


(2.5) 


— OO 

OO 


/ 


/ OO \ 

E^{^l\u)+ [ Ef(A|n)H(A,/i)dA 


\E 


R 


(A))o(i + h) 


E^{^l\u)+ I E^{X\u)V{X,^l)dX 


R( 


( 2 . 6 ) 


/ 
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Here I is the identity operator in Ti. The action on the bra-vectors {E^{X)\ is quite similar. Below 
we shall mark the operators of this type by “tilde” and denote their action by the sign “o”. 

The second type of operator acts in the space fi. It is convenient to treat these operators as 
2x2 matrices with operator-valued entries: 


/ 


A = 


An(u,v) Ai 2 (u,v) 


y A2i(u,v) A22(u,v) 


The action of the operators (2.7) is given by 


/ 2 °° \ 
^ [ Aik{u,v)Ej^{X\v)dv 


A-\E^{X)) = 


k=l_ 


/ A 2 k{u,v)Ej^{X\v)dv 


fc=i_ 


/ 


(2.7) 


( 2 . 8 ) 


(£’'^(A)| • H = j ^ f Ef{X\u)Aji{u,v)du;Y^ f Ej{X\u)Aj 2 {u,v)du\ (2.9) 

i=i-bo / 

This type of operator will be marked by “hat”. 

We would like to draw to the attention of the reader that we use “hat” not only for the operator¬ 


valued matrices of type (|2.7D , but for their matrix elements also. We hope, that this system of 
notation will not cause misunderstanding. 

Finally, we define the trace of “hat”-operators as 


TrA = 


OO 

tr Hii-I-tr H 22 = J (^Aii{u,u) + A 22 iu,u)^ du, 


( 2 . 10 ) 


where 


tr Ajk = 


Ajk{ 


u, u 


du. 


( 2 . 11 ) 


2.2 Dual fields 

The important objects, on which the r.h.s. of ( ^.1| ) depends, are auxiliary quantum operators — dual 
fields, acting in an auxiliary Fock space. They are introduced in order to remove two-body scattering 
and to reduce the model to free fermionic. One can find the detailed definition and properties of 
dual fields in Section 5 and Appendix C of |^. Here we repeat them in brief. 
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Consider an auxiliary Fock space having vacuum vector |0) and dual vector (0|. The three dual 
fields tpiX), (pDiiX) and (pA 2 {X) acting in this space are defined as 


4>A2{X) = qA2 {X)+pd2{X), 

</>Di(A) = gDi(A)+ pai(A), (2.12) 

V'(A) = q^p{X)+Pij{X)- 

Here p{X) are annihilation parts of dual fields: p(A)|0) = 0; q{X) are creation parts of dual fields: 
(0|g(A) = 0. Thus, any dual field is the sum of annihilation and creation parts (the dual field ip{X) 
should not be confused with the field which appears in the expression for the Hamiltonian 

of the model). 

The only nonzero commutation relations are 


\PAi{X),q^{p)] = \Pi,{X),qA 2 {^^)] =ln/i(/x,A), 
\PD 2 {X),qij{p)] = \p^{X),qDA^^)] =ln/i(A,/i), 
\Pi^iX),q-,pip)] =ln[h{X,p,)h{p,X)], where h{X, p) = 


(2.13) 


X — p + ic 


ic 


Recall that c is the coupling constant in (1.2). It follows immediately from (2.13) that the dual fields 
belong to an Abelian sub-algebra 


bP{X),'lp{p)] = [lp{X),cl)a{p)] = [0fe(A),())a(^)] = 0, 


(2.14) 


where a,b = A 2 ,Di. The properties (p.l3 ), ( 2.14 ), in fact, permit us to treat the dual fields as 
complex functions, which are holomorphic in some neighborhood of the real axis. 


2.3 Fredholm determinant 

Now we are ready to describe the determinant representation (^.Ij). 


The important factor entering the r.h.s. of (^) is the Fredholm determinant of the integral 
operator I -\-V, acting on the real axis as 


(i+ V) o f{p) = f{X) + f V{X,p)f{p)dp, 
^ ' J—oc 


(2.15) 


where /(A) is some trial function. Recall that I is the identity operator in the space Ti. 
The kernel V (A, p) can be presented in the form 

(E^(X)lF^(p)} 


V(X,p} = 


X — p 


(2.16) 
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where the components of vectors (-E^(A)| and |£'^(/u)) are equal to 

E^{X\u) = E^{X\u) = E+{X\u), 
E^{X\u) = -Ef{X\u) = E_{X\u). 

Due to this specification we have 


(2.17) 

(2.18) 


(i^^(A)|i7^(A)) = 0, (2.19) 

hence the kernel V{X,fi) possesses no singularity at A = /r. Later on, the orthogonality property 


( 2.1£ ) will play an important role. 

The functions E± introduced in Q are equal to: 


^ , 1 Z(u,X) / e-^^2 

E+{X\u) = —I-+ 




27r Z{u,u) \u — X + iO u — X — iO 


m 


xe 


V'(M) + T(M)+i((/>i3, (\)+(l>An (A)--i/>(A)-r(A)) 


( 2 . 20 ) 


L;_(A|u) = ^Z(u,A)e5('^^iW+'^^2W-V’W-r(A)) 

ZTT '' 

where the function Z{X,fi) is defined by 

Z(A,/i) = 


g-<i>Oi(A) g-0/i2(A) 


h{fi,X) h(A,/i) 
Here ip{X), 4 >DiW and 4>a2W are just the dual fields ( p.l2 ). 
Recall also that the function ??(A) is the Fermi weight 

re(A)--' 


'd(A) = I 1 + exp 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


This function defines the dependence of the correlation function on temperature and chemical po¬ 
tential. The function t(A) is the only function depending on time and distance: 

r(A) = RA^ — ixX. (2.24) 

Thus, functions E±{X\u) depend on time t, distance x, temperature T and chemical potential h, but 
this dependence, as a rule, suppressed in our notations. 

In order to define the function bi 2 {u,v), we introduce vectors (F^(A)| and \F^{fi)), belonging to 
the space 7i 

' Fi^(/x|u) ^ 

(F^(A)| = (t/'(A|u),F2^(A|u)); |F^(f^)) = , (2.25) 

V T2^(/x|u) y 
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as solutions of the integral equations (see [Q]) 


OO 

(/ + y)o(F^(^)|^(F^(A)|+ J V{X,^^){F^{^^)\d^={E^{X)\, 

— OO 
OO 

|F^(A))o(/ + f) = |F^(^))+ J \F^{X))V{X,fi)dX = \E^{fi)) 

— OO 

Let us define operator B as 

OO OO 

B= j \E^{X)){E^{X)\dX= J \E^{fi)){E^{fi)\dfi. 

— OO —OO 

The components of this operator are 

OO 

Bjkiu,v)= J T/^(A|u)L;fc (A|u)dA, j,k = 1,2. 

— OO 

The function 612 (u,u) is the matrix element of the operator b defined by 

b = B + g. (2.30) 

Here g is 

gi2iu,v), 5i2(n,u) = (5(tt - (2.31) 

have 

= Bab{u, v) for all a, b except a = 1, 6 = 2, (2.32) 

= Bi2{u,v) - gi2{u,v) (2.33) 

The importance of operators b and g will be clear later, after the formulation of the Riemann-Hilbert 
problem. 

The last factor in the representation of the correlation function, which is not explained yet, is 
the determinant of the operator I — ^Kt- This operator also acts on the whole real axis. Its kernel 
is given by 

grO.f.) = ( ^, ^^, ) vTaIvTaO. (2.34) 

This determinant does not depend on time and distance, therefore it can be considered as a constant 
factor. 



Due to the above definition we 

bab{u,v) 

bi 2 {u,v) 


(2.26) 

(2.27) 

(2.28) 

(2.29) 





Thus, we have described the r.h.s. of (2.1). The temperature correlation function of local 
fields is proportional to the vacuum expectation value in the auxiliary Fock space of the Fredholm 
determinant of the integral operator. The auxiliary quantum operators—dual fields—enter the 


kernels V and bi 2 {u,v). However, due to the property (2.14) the Fredholm determinant is well 
defined. 


3 Riemann—Hilbert problem 

Consider the Riemann-Hilbert problem for the operator x('^) acting in the space Ti. and depending 
on the complex parameter A; 


x(A|u,u) = 


Xii(A|u,u) xi2iMu,v) 

X2 i(A|u,u) X22(A|u,u) 


(3.1) 


The matrix x{X\u,v) is holomorphic for A G C \ R, and satisfies the normalization condition at 
A = oo 


, 10 , 

x{oo\u,v) = I = (5(u - u). 

0 1 

The boundary values on the real axis are related by 

X-(A) = x+(A)G(A), AgR, 

where we set X±(A) = lime^+ox(A ± ie)- The jump matrix G'(A) here is given by 

G{X\u,v) = ( \ 6{u — v) + 27riH{X\u,v), 


where 


HiX\u,v) = \E^{X)){E^{X)\ = 


' Hf(A|u)F;f(A|u) E{^{X\u)E^{X\v) 
y E^{X\u)E^{X\v) E^{X\u)E^{X\v) 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


Recall that due to (RlTf) , (12.Iq ) 

(A|u) = ii;+(A|u), E^{X\u) = E.{X\u), 

E^{X\v) = -E_iX\v), E^iX\v) = E+iX\v), 

where the functions Fi±(A|M) are defined in (2.2C) and (2.21). However, we do not use the explicit 


(3.6) 


expressions for E± in this section. 
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One should understand the r.h.s of the jump condition 
TL. Another words, more detailed the equality (3.3) means 

2 

^~^—oo 

Thus, we are dealing with an infinite dimensional Riemann-Hilbert problem. 
The orthogonality property ( p.l9 ) 

{E^{X)\E^{X)) = 0 

implies TrH^{X) = 0 for n > 1 . Hence, 

detG(A) = 1. 


as an operator product in the space 


(3.7) 


(3.8) 


(3.9) 


We assume that the Riemann-Hilbert problem is solvable. We note also that in the class of 
integral operators which we are dealing with, the analyticity of the operator with respect to the 
parameter A implies the analyticity of its determinant. Therefore, due to the Liouville theorem and 
), we conclude that detx(A) = 1. Hence, as in to the usual matrix case the 


equations (3.2), 

solvability of our Riemann-Hilbert problem implies the uniqueness of the solution. It means also 
that the inverse matrix x“^(A) exists and is holomorphic for A G C \ R. 

The integral equations (|2.26|) , (^.271) 


{l + v)o{E^{^^)\ = {E^{X)\, 
\F^iX))o[i + v)=\E^ifi)), 


(3.10) 

(3.11) 


were used in |0] for derivation of the differential equation, describing the Fredholm determinant. 


Now we prove the equivalence of these integral equations and the Riemann-Hilbert problem (3.2), 
(|3.3| ). Namely, we establish that the solution of the Riemann-Hilbert problem can be presented in 
terms of the solutions of integral equations and vice versa. 

The following integral representation is the basis of this paper. 


Theorem 3.1 A solution of the Riemann-Hilbert problem (3.5) has the following integral repre¬ 
sentation 


m = i- 


\E^{pi)){E^{p.)\ 

pi — X 


dh, 


(3.12) 


where \E^{pL)) is the solution of the integral equation (3.11): 


?R 


{X))o[l + v)=\E^{fr)). 
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The inverse of this solution has the following integral representation: 


r'(A) = / + 


fj, — X 




(3.13) 


where (F^(A)| is the solution of integral equation MM)- 

(i+v)o{F^{g:)\ = {E^{X)\. 


Proof. The operator x(A) defined in ( |3.12| ) is a holomorphic function everywhere except the real 
axis and possesses the correct asymptotic behavior when A —> oo. On the other hand, on the real 
axis we have 


X+(A)G(A) = x+(A) + 27 ri 


E^{X)){E^{X)\ 


\F^{^i)){E^{^i)\E^{X)){E^{X)\ 


fj, — X — iO 


dp, 


= X+(A) + 27ri 


|S«(A))- / \F^{p))V{p,X)dp 


R/ 


{E\X)\ 


= I - 


= I - 


\F^{p)){EL{p)\ 
p — X — if) 

\F\p)){E^{p)\ 
p — X + if) 


dp + 2T:i\F^{X)){E^{X)\ 


dp = X-(A). 


(3.14) 


Thus, the operator-valued matrix x(A) satisfies all the conditions of the Riemann-Hilbert problem 

g), (P). 

Just the same method allows one to check representation (|3.13[). The theorem is proved. 


A direct corollary of the representation (3.12) is the following asymptotic expansion of X±(A) at 
A —> oo ^ ^ 

X±(A) = I + — + ^ + ..., (3.15) 

where 

OO OO 

B= f \F^{X)){E^{X)\dX, C= f X\F^{X)){E^{X)\dX, ... (3.16) 


We see that the operator B introduced in the previous section (see (p.28)) appears as the first 


coefficient in the asymptotic expansion (3.15) 
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Theorem 3.2 For arbitrary A € C 


(3.17) 

(3.18) 


|F«(A)) = x(A)|i5;^(A)), 

(F^(A)I = (ii;^(A)|r'(A). 

In particular, for A G R 

|F«(A)) = x+(A)|i7^(A)) = x-(A)|ii^^(A)), (3.19) 

(F^(A)I = (ii;^(A)|xT'(A) = (i?^(A)|x:'(A). (3.20) 

Proof. Using the definition of vectors |F^(A)) and (F'^(A)| and integral representations for x(A) 
and X~^(A) we have 

—oo 
oo 

= |ii;''(A))- J \F^{^i))V{fi,X)dp=\F^iX)), (3.21) 

— OO 

—oo 
oo 

= {E^{X)\- J V{X,p){F^{^r)\d^r = {F^{X)\. (3.22) 

— OO 

Due to the property (ii^^(A)|U'^(A)) = 0, we find for A G R 

X-(A)|ii;^(A)) = x+(A)(/ + 27 ri|ii;«(A))(i?^(A)|)|ii;^(A)) 

= x+(A)|i?^(A)), (3.23) 

(U^(A)lxl'(A) = (ii;^(A)|(/-27ri|ii;^(A))(ii;^(A)|)x;'(A) 

= (ii;^(A)|x;'(A). (3.24) 

Here we have used 

G-i(A) = / - 27 ri|.U^(A))(.U^(A)|. 

The theorem is proved. 

The latest theorem allows one to consider the transformation 

|i^«(A)) ^ |F«(A)) = x(A)|i7^(A)) (3.25) 

as a gauge transformation, which is analytic in the complex plane A G C. 
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Thus, we have proved the equivalence of the integral equations ( |3.1CI| ), (|3.11 ) and the Riemann- 
Hilbert problem ( |3.2|) , . It is interesting to mention that integral equations allows us to consider 

the transformation |T’'^(A)) —> |i?^(A)) as transformation in the space TC. At the same time, the 
solution x(^) of the Riemann-Hilbert problem defines the same transformation in the space Tl. 


4 Lax representation 

In the previous paper Q we had obtained a generalization of the Nonlinear Schrodinger equation, 
which describes the Fredholm determinant det (i + V) 


-idtbi2{u,v) = -dlhi2{u,v) 

OO 

+2 J dwidW2bi2iu,Wi)b2liwi,W2)bi2{'W2,v), 

— OO 

idtb2iiu,v) = -dlb2i{u,v) 

OO 

+2 / dwidW2b2liu,Wi)bi2iwi,W2)b2l{'W2,v). 


(4.1) 


(4.2) 


In this section we derive the equations (0), (PD and their Lax representation in terms of the 
Riemann-Hilbert problem (|3.3[ ). Namely, we prove that Riemann-Hilbert problem (^]^) implies 
equations (fl.lD , (P)- 

In the paper |Q] we have obtained the system of linear differential equations for vectors 
and (£^^(A)|. Let us recall this system 


= L{f,)\E^{fi)), a,(S^(A)| = -(F;^(A)|L(A), 

dtlE^ifi)) = M(m)|L;^(m)), dt{E^{X)\ = -{E\X)\MiX), 


where we have set 


L(A) — Aft + [ 5 , fj], M{X) ——XL{X) + dxg- 


(4.3) 

(4.4) 


(4.5) 


Here the operator g is defined in (2.31), the operator d is given by 

. 1/1 0 


^ 2i 


0 -1 


5{u — v). 


(4.6) 
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The proof of the system (^^), (^~^) is based on the obvious relations 

(4.7) 

It is easy to check that the system (|4.3|) , (|4.4|) is compatible. Indeed, the compatibility condition 

dti - d^M + [L, M] = 0 


implies 


which in turn is equivalent to 


[dtg, a] - dig + 


[5, d] ■ d^g 


= 0 , 


(4.8) 


idtgi2iu,v) - dlgi2{u,v) = 0 . 


The last equality is obviously valid due to (|2.3l[) . 

Our aim now is to find the system describing the derivatives of vectors |T^(A)) and (T^(A)| with 
respect to x and t. The solution of this problem is given by the following theorem. 


Theorem 4.1 The vectors |F^(A)) and {F^{\)\ satisfy the following system of differential 
equations 


9,|T^(A))=T(A)|T^(A)), 

d,(F^(X)l = -(F^(X)l£(X), 

(4.9) 

9i|F^(A)) = MA)|F^(A)), 

dt(F^(X)j = -(F^(X)lM(X), 

(4.10) 

where we have set 



C{X) = 

Aft + [b,a], 

(4.11) 

M(X) = 

—XC{X) + dxb. 

(4.12) 


Here the operators b and a are defined in ( 2. 3(\) and (4-6). 

Proof. It follows immediately from ( p.l7|) , (|3.1^ ) and (4.3), (4.4) that the derivatives of vector 
|F^(A)) with respect to x and t can be presented in the form 

a,|F«(A)) = T(A)|F«(A)), 
dt\F^{X)) = M{X)\F^{X)), 


(4.13) 

(4.14) 


where 

C = d^xW ■ x~^{X) + x(A) • L(A) • 

M = dtxW • r'(A) + X(A) • M(A) • x-'(A). (4-15) 
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It was shown in the previous section that the gauge transformation ( 3.25| ) is continuous across the 
real axis. The direct corollary of this property is that matrices L and Ad possess no cuts on the real 
axis. Let us prove this directly. 

Let = lime_>+oT(A ± ie). Then we have for ImA = 0 


^±(A) = a,x±(A) • xi'(A) + X±(A) • L(A) • xi'(A). 

Using the jump condition (3.3) we have 

T_(A) = a,x-(A)-xl'(A)+X-(A)-L(A)-x:'(A) 


(4.16) 


= a,x+(A) • x+ (A) + x+(A) • d^G{\) • G-^(A) • x+ (A) 
+ x+(A)-G(A).L(A)-G-1(A)-x;'(A). 

Via (1^) , (^ ) and due to (^(4|) , (|3(5|) we find 


a,G(A) = [L(A),G(A)]. 


After the substitution of this expression into (1.17) we arrive at 

T_(A) = 9,x+(A) ■ x;'(A) + x+(A) • L(A) • x;'(A) = 4(A). 


(4.17) 


(4.18) 


Thus, L is continuous across the real axis. One can perform just the same consideration for matrix 
M. also. 

Due to the analyticity properties of x(A), X~^(A) and their normalization conditions, we conclude 
that matrices £(A) and Ad (A) are holomorphic for A G C and possess the asymptotics 


T(A) ^ A<t + 0(1), M(A) ^-A2d + C>(A). 


(4.19) 


Due to the Liouville theorem we conclude that >C(A) is a linear function of A, while Ad (A) is a 
quadratic function of A. Using the asymptotic expansion ([1.15 ) we find 


£(A) = Ad + [6, d], 

Ad (A) = — A^d — \[b, d] + dxQ + [d, C] + [b, a] ■ B — B ■ [g, d]. 


(4.20) 

(4.21) 


The formula ( 4.20D exactly reproduces the expression (4.11). In order to reduce the expression for 
Ad (A) to the formula (4.12), one should take into account that relations (|4.15|) provide us an infinite 
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set of identities between decomposition coefficients B, C,... and their derivatives with respect to x 
and t. For exampfe, the first of the refations (4.15) implies 


d.x{\) = C{X)x{\) - mm- (4.22) 

One can substitute the asymptotic expansion for x(^) and explicit expressions for C{X) and L{X) 
into the last formula. After this, comparing coefficients at negative powers of A, we obtain the 
mentioned above set of identities. In particular, for we have 

d,B = [a, C] + [b,a]-B-B-[g,a], (4.23) 


Therefore, we arrive at 

M{X) =-X‘^a-X[b,a]+d^b. (4.24) 

Thus, we have proved differential equations (|4.9|) , ([4.1CI|) for the vector |F^(A)). The equations for 
vector (T^(A)I can be found in a quite similar way. 

Since a gauge transformation does not disturb the compatibility condition, we obtain for the pair 
and ( |4.10| ) 

dtC - dm + [A M] = 0, 


what implies in turn 


[dtb, a] - dlb + 


[b, a] ■ d^b 


= 0 . 


(4.25) 


The matrix operator-valued equation (4.25) is equivalent to the four scalar operator-valued partial 
differential equations. It is easy to check that the equations for the diagonal part of b are valid 
automatically due to identity (|4.25): 


dxbii{u,v)=i / bi2iu,w)b2i{w,v)dw, 


dxb 22 iu,v) =-i / b 2 i{u,w)bi 2 {w,v) dw 


(4.26) 


These equations, being substituted into the antidiagonal part of (4.25), give the non-Abelian Non¬ 
linear Schrodinger equation ([4.1D, ([4.2|). More compactly they can be written in the following form 


- idtbi 2 = — d ^ hi 2 + 2612621612, 


idthi = -< 9^621 -I- 2 A 1612621 . 


(4.27) 
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One should understand here the nonlinear terms in the sense of integral operator products. 

Remark. The method of derivation of the Nonlinear Schrodinger equation described in this 
section is closely connected with the “dressing procedure” proposed in 


5 Modification of the Riemann—Hilbert problem 


In the previous sections we have demonstrated that integral operator I + V generates in a natural 
way the Riemann-Hilbert problem. The latter, in turn, defines the dressing gauge transformation, 
which permits us to obtain the exactly solvable classical differential equations. As was shown in the 
paper Q, the Fredholm determinant of the operator I + V can be described in terms of the solutions 
of the Riemann-Hilbert problem and differential equations mentioned. Namely, the logarithmic 
derivatives of the determinant with respect to distance and time are expressed in terms of the 
operator x asymptotic expansion coefficients B and C. Let us present here the list of logarithmic 
derivatives (see @): 

dx log det+ = ztrRii, 

dtlogdet(i + V) =itr{C22 - Cii - ^21512), 

. .. _ ( 5 - 1 ) 

5^5a;logdet(/ + R) =-tr(Si2.B2i), 

dtdx logdet^I + r) = itr{dxbi 2 ■ B 21 - dxB 2 i ■ bu). 

Thus, the solution x of the Riemann-Hilbert problem allows to reconstruct the Fredholm deter¬ 
minant up to a constant factor, which does not depend on x and t. Hence, the calculation of the 
correlation function of local fields is reduced to the solving of the operator-valued Riemann-Hilbert 
problem. 

In the present section we formulate a new Riemann-Hilbert problem, which appears to be more 
convenient from the point of vie of an asymptotic analysis. The detailed asymptotic investigation of 
this problem will be given in our forthcoming publication. Here we restrict our selves only to basic 
formulations. 

Consider the following substitution 

x(A) = x(-)(A)/(A), (5.2) 

where the triangular matrix y;^(A|rt, u) is defined by 

n ( 1 0 \ / 0 -1 

X (A|u,u) = h(u-u) -h 

01/ \ 0 0 


912[U,V) 

u — X 


(5.3) 
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The expansion of at A —*• oo is given by 


F(A)=/-^- 


9 idxQ 


+ 


A A2 

Due to ( 3.15 ) we find for the asymptotic expansion of 

x^^\x) = i + \ + ^ + ..., 

where the operator-valued matrix b was defined in ( 2.301 ): b = B + g; and 


(5.4) 


(5.5) 


c = C + Bg + idxg. 


(5.6) 


The new operator-valued matrix x^”^^ (A) satisfies the modified Riemann-Hilbert problem: 

A^CX), 


(5.7) 

£™)(A) = xf^(A)G(™)(A), AgR. 

The modified jump matrix is equal to 

GM(a)=x°(A)G(A)(x°)-'(A), (5.8) 

and possesses the following entries 

G^^\x\u,v) = 6{u-v)-6{u-X)Z{v,X)'d{X)e‘^^i^^\ (5.9) 

G^^\x\u,v) = -27rz(l - 79(A))<5(n - A)5(n - A)e^(^)+^(^) (5.10) 

G^^\x\u,v) = r9(A)Z(M,A)Z(n,A)e'^^i(^)+‘^^2W-V’D)-RA)^ 

27r 

G^^\x\u,v) = 6{u-v) - Z{u,X)6{v - X)^{X)e'f‘^2W_ (5.12) 


To derive the above matrix elements, we have used the relation 

Z{u, u) = -h e "'^^2 (“). 

Remark. One should understand the equations (^)^ ) - (|5.12 ) 
Namely, 


in the weak topology sense. 


Hi{ui,U2) = H2 {ui,U2), 

OO 

j fl{ui)Hi{ui,U2)f2{u2)dUidU2 

— OO 

where fi{ui),f 2 {u 2 ) are test functions. 


{Weak topology) 

OO 

/ fl{ui)H2iui,U2)f2iu2)dUidU2. 
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The transformation, considered above, is the direct generalization of the approach proposed in 
1^] for the free fermionic limit of the quantum Nonlinear Schrodinger equation. The main advantage 
of the modified Riemann-Hilbert problem ( ^.T] ) is the simple explicit dependency of the jump matrix 
on the variables x and t. This allows us, in particular, to use another method for deriving 
the differential equations (4.27). Indeed, the jump condition for the matrix 


X^’^\X\u,v) = x(”")(A|'u,u)exp ('-^r(A)cJ3 


can be written as 


£^)(a) = xJ)(a)g(^)(a), 


(e)/ 


A G R, 


where the jump matrix 

G(®)(A) = exp (^-^r(A)cJ 3 ^ G(”*)(A)exp Qr(A)(T 3 
does not depend on x and t. Hence, the logarithmic derivatives 

-^x(A) = (a,x(^)(A)) (x(^))-'(A), 


(5.13) 

(5.14) 

(5.15) 

(5.16) 

(5.17) 


possess no cut on the real axis. Therefore they are holomorphic for A G C and have the following 
asymptotics: 

^ Ad, J^t{X) ^ -A^d. (5.18) 


Due to the Liouville theorem we conclude that J-x{X) and J~t{X) are linear and quadratic functions 
of A respectively. Using the asymptotic expansion (^]^) we arrive at 

Tx{X) = Ad+[S,d], (5.19) 

jr^(A) = -A^d-A[S,d]+ 9,,6. (5.20) 


Thus, we again have obtained the Lax representation ( [4.11| ), (|4.12 ). 

Finally, let us rewrite the logarithmic derivatives of the Fredholm determinant in terms on new 
decomposition coefficients h and c. We have 


Sa; log det^I + =itrSii, 

dxdxlogdet(^I + V^ = -tr (612621), 

logdet(^/ +='itr(c22 - cii), 

9 f 9 a,logdet(/ + =itT{dxbi2 ■ 621 - dxhi • ^12)- 


(5.21) 
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Summary 


The main purpose of this paper was the formulation of the Riemann-Hilbert problem associated with 
the correlation function of the quantum Nonlinear Schrodinger equation. We used the Riemann- 
Hilbert problem in order to derive the non-Abelian classical Nonlinear Schrodinger equation, which 
describes the Fredholm determinant. As we have seen, the solution of this equation is completely 
described by the solution of the Riemann-Hilbert problem. This permits us to reduce the calculation 
of the Fredholm determinant to the investigation of the Riemann-Hilbert problem. The detailed 
asymptotic analysis of the latter will be performed in our forthcoming publication. 
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